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Abstract. Cravitational theories with multiple scalar fields coupled to the metric and each 
other — a natural extension of the well studied single-scalar-tensor theories — are interesting 
phenomenological frameworks to describe deviations from general relativity in the strong-field regime. 
In these theories, the A-tuple of scalar fields takes values in a coordinate patch of an A-dimensional 
Riemannian target-space manifold whose properties are poorly constrained by weak-field observations. 
Here we introduce for simplicity a non-trivial model with two scalar fields and a maximally symmetric 
target-space manifold. Within this model we present a preliminary investigation of spontaneous 
scalarization for relativistic, perfect fluid stellar models in spherical symmetry. We find that the 
scalarization threshold is determined by the eigenvalues of a symmetric scalar-matter coupling matrix, 
and that the properties of strongly scalarized stellar configurations additionally depend on the target- 
space curvature radius. In preparation for numerical relativity simulations, we also write down the 
3 + 1 decomposition of the field equations for generic tensor-multi-scalar theories. 
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1. Introduction 


Modifications of general relativity (GR) often lead to the introduction of additional 
degrees of freedom . The simplest and best studied extension of GR is scalar-tensor 
(ST) theory, in which one or more scalar fields are included in the gravitational sector 
of the action through a non-minimal coupling between the Ricci scalar and a function 
of the scalar field(s). A further motivation to study ST theories is that they appear in 
different contexts in high-energy physics: they can be obtained as the low-energy limit 
of string theories , in Kaluza-Klein-like models or in braneworld scenarios . 
Moreover, ST theories play an important role in cosmology [^. 

Almost 60 years ago, in an attempt to implement Mach’s ideas in a relativistic 
theory of gravity, Jordan, Fierz, Brans and Dicke proposed a specific ST theory 
(commonly referred to as “Brans-Dicke theory”) as a possible modification of GR [7||^ . 
Their theory is still viable, but it has since been constrained to be extremely close 
to GR by Solar System and binary pulsar observations [^. Brans-Dicke theory 
was generalized by Bergmann and Wagoner, who considered the most general ST 
theory with a single scalar field and an action at most quadratic in derivatives of the 
fields innii- In 1992 Damour and Esposito-Farese introduced and investigated tensor- 
multi-scalar (TMS) theory, a generalization of the Bergmann-Wagoner theory to an 
arbitrary number of scalar fields |13| . Multiple scalar degrees of freedom are a generic 
prediction of string theories and theories involving extra dimensions (see e.g. |14[|17| ). 
In recent years, Bergmann-Wagoner theory has been extensively studied in the case of a 
single scalar field (see e.g. 18 -21 and references therein). In comparison, very limited 
attention has been devoted to the phenomenological implications of TMS theory. 

Even in the simplest case of a single scalar, ST theories give rise to interesting 
phenomenology. Although their action is linear in the curvature tensor, and scalar- 
matter couplings are highly constrained by observational bounds from the Solar 
System j^, ST theories can modify the strong-field regime of GR. Indeed, the 
equations of structure describing compact stars can admit non-perturbative solutions 
where the scalar fields can have large amplitudes [^. This phenomenon, called 
spontaneous scalarization, can significantly affect the masses and radii of neutron stars. 
Spontaneous scalarization is strongly constrained by binary pulsar observations |24| , 
but it could still leave a signature in the late inspiral of compact binaries through 
the so-called “dynamical scalarization” 25 - 27 , which may be observable by advanced 


gravitational-wave detectors 29 

In this article we study the phenomenology of TMS theories. The importance 
of our work for black hole physics is to a significant extent of indirect nature. 
Stellar collapse represents, of course, one of the most important channels for black 
hole formation, and ST theories of gravity are more likely to produce experimental 
signatures during black hole formation than in the dynamics of the remnant “quiescent” 
black hole spacetimes 30 31 . The reason is that there are strong no-hair theorem^ 


implying that stationary, vacuum black hole solutions in ST and TMS theories must be 
identical to GR (cf. [O 34-38 and Sotiriou’s contribution to the present volume (^). 


In addition, it has been shown that the dynamics of black-hole binaries in ST theory 
is undistinguishable from that in GR up to 2.5 post-Newtonian order 40,41 or - in 
the case of extreme mass-ratio - to all post-Newtonian orders [42]. This result does 


* Hairy black hole solutions are possible in presence of a complex scalar field (i.e. two real scalar fields), 
as long as its phase is time-dependent [32] . No-hair theorems can also be evaded if the potential is 
non-convex; in this case the solutions are unstable, but their growth time can be extremely large [33] . 
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not apply in the presence of non-trivial boundary conditions 43 ■ 45 , but in this case 
violations of the no-hair theorem would most likely be unobservable. Nevertheless, 
black hole binaries have been studied in the framework of single-scalar theories using 
numerical relativity (cf. |44{[45] ), and the extension of such studies to TMS theories 
may still reveal surprising discoveries. 

In order to pave the way for numerical investigations of black holes and neutron 
stars in TMS theories, in Section]^ we write down the TMS field equations in a 3+1 
formalism that is suitable for numerical evolutions of compact binary systems. 

A promising avenue to understand the experimental implications of TMS theories 
(following the reasoning of (23p^ ) is to focus first on the coupling between the various 
scalars and matter. The scalar fields take values in a coordinate patch of a Riemannian 
target-space manifold. It is natural to ask whether this additional geometric structure 
can leave a detectable signature in compact stars and/or in the late inspiral of compact 
binaries, while still allowing for solutions compatible with binary pulsar observations. 

This paper is a preliminary investigation in this direction. We will mainly 
focus on a simple non-trivial TMS model with two scalar fields and a vanishing 
potential. This model contains the main novel features that distinguish TMS theories 
from single-scalar theories, i.e. the presence of a target-space manifold with non¬ 
trivial Riemannian structure and non-vanishing curvature, and the presence of a 
continuous symmetry that is spontaneously broken by scalarized solutions. A more 
systematic study of scalarization in this prototype TMS theory will be the topic of 
a future publication, and it should give us useful insight into strong-field effects that 
characterize more general TMS theories. 

The plan of the paper is as follows. In Section we introduce the action and 
the field equations of TMS theory. We then specialize the field equations to the case 
of two scalar fields with vanishing potential and a maximally symmetric target-space 
manifold. In Section we derive the equations for slowly rotating relativistic stars, 
and we perform numerical integrations of these equations in the non-rotating case. In 
Section we derive the 3 + 1 decomposition of the TMS field equations, which can 
be used to perform fully non-linear numerical evolutions. In Section we draw some 
conclusions and point out possible directions for future work. The Appendices contain 
some technical material on the structure of target spaces in our model ( [Appendix A ), 
experimental constraints (Appendix Bl and perturbative arguments to predict the 


spontaneous scalarization threshold in models with two scalar fields (Appendix C \. 


2. Tensor-multi-scalar theories: action, field equations, scalar-matter 
couplings, and symmetries 

2.1. Units, notation and conventions 

Throughout the paper we use units with c = 1. The gravitational constant measured 
in a Cavendish experiment is denoted by G, while the “bare” gravitational constant 
appearing in the action is denoted by G*: the relation between the two constants 


is written down explicitly in Eq. (B.l). Indices on space-time tensors are denoted 


by Greek letters and take values 0,..., 3, and space-time coordinates are denoted by 
x^. The Lorentzian space-time metric is taken to have signature (—,+,+,+), and 
its components are denoted by The conventions for the Riemann curvature 

tensor and its contractions, as well as the notation for symmetrization and anti- 
symmetrization of tensors, are those of Misner, Thorne, and Wheeler (T^. 
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The A^-tuple of scalar fields ..., takes values in a 

coordinate patch of an iV-dimensional Riemannian target-space manifold. Indices 
on target-space tensors are denoted by early capital Roman letters A, B, C, .. and 
take integer values 1,..., fV. Components of the Riemannian target-space metric are 
denoted by "fAsisp), and the associated Christoffel symbols are denoted by 'y^ABif)- 
The target-space Riemann curvature tensor is denoted by with obvious 

notation for derived quantities such as the Ricci tensor and the Ricci scalar. If the 
target space has a Hermitian structur^ then indices on complexified tensors are 
denoted by lower-case Roman letters, and take values 1,2,..., N/2. Holomorphic 
coordinates are denoted by and the components of the Hermitian metric in 

these coordinates are denoted by 'jabift'f)- For reference, in Table we provide an 
overview of the meaning of the various symbols and conventions used in this paper. 


2.2. Action and field equations for N real sealars 

We consider a gravitational theory with metric tensor and scalar fields tp^,..., (p^ 
which take values in a coordinate patch of an TV-dimensional target-space manifold. 
We assume that all non-gravitational fields, denoted collectively by dt, couple only to 
the Jordan-frame metric = A^{{p)g^,j, so that the matter action has the functional 
form This assumption guarantees that the Weak Equivalence Principle 

(WEP), which has been experimentally verified with great accuracy [^, will hold. 
The quantity A(ip) is a conformal factor relating the metrics and g^,^. 

The most general action which is invariant under space-time and target-space 
diffeomorphisms (up to boundary terms and field redefinitions), and has at most two 
space-time derivatives, can be written in the form |13| 


S = 



^g^''-fAB{‘p)di,(p^ d^(p^ 


Viq>) 




( 2 . 1 ) 


where G* is a bare gravitational constant, and g and R are the determinant and Ricci 
scalar of respectively. The positive-definiteness of the target-space Riemannian 
metric ‘jABi^p) guarantees the absence of negative-energy excitations. The scalars 
are dimensionless and the potential V{(p) has length dimensions minus two. The 
conformal factor A{ip) is dimensionless. In the case of a single scalar (iV = 1), the 
target-space metric jABip) reduces to a scalar function jip), and the choice ■j^p) = 1 
can be made without loss of generality. 


The field equations of the theory, obtained by varying the action (2.1 1 with respect 
to g^" and ip, take the form 

= 2jab{p)'^+ ‘^V{p)g^i, + SttG* - ]^Tg^^ , (2.2) 

Uip^ = - 47rG„7^^ {p) ^ T. 


dpi 


dpi 


(2.3) 


Here is the covariant derivative associated with g^^, and □ = is the 

corresponding d’Alembertian operator. The Ricci tensor built out of the metric ( 7 ^^, 


* For an introduction to Hermitian structures and complex differential geometry, see e.g. |48 
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G 

G* 

Vm 

A,B,C 

N 

lAB 

Ibc 

R^bcd 

4- 

A{ip) 

V{ip) 
a, &, c 

^ah 

hc^ /6c 

r 

* — * 
a , a 

6 

a,/3i 

Z 

w 

Gravitational constant from a Cavendish experiment 

Bare gravitational constants appearing in the action 
Spacetime indices 

Spacetime coordinates 

Spacetime metric in the Einstein frame 

Covariant derivative associated with 

Spacetime Riemann tensor 

Scalar-field indices in real notation 

Number of scalar fields 

Gravitational scalar fields in real notation 

Target-space metric in real notation 

Christoffel symbols on the target space in real notation 
Target-space Riemann tensor 

Non-gravitational fields 

Spacetime metric in the Jordan frame 

Einstein-Jordan frame conformal factor 

Scalar-field potential 

Scalar-field indices in complex notation 

Gravitational scalar fields in complex notation 
Target-space metric in complex notation 

Christoffel symbols on the target space in complex notation 
Target-space curvature radius (N = 2) 

Scalar-matter coupling function (TV = 2), see Eq. (2.191 
Linear-term coefficients in the expansion of log A^ipj’ijj) 
Quadratic-term coefficients in the expansion of log A{ij}, pj) 
Generic rotation angle in the target-space complex plane 
Redefinition of a*, after rotation 

Redefinition of the fields 75 , (p after rotation 

Re [■)/'] 

Im ['(/;] 

t,r,e,(j) 

v{r) 

m(r) = rfi(r) 
uj{r) 

n 

p 

p 

Ub 

Subscript “0” 
R, R 

M 

Q 

Mb 

K, no, mb, 7 

Spacetime coordinates for stellar models 

Lapse function 

Mass function 

Fluid differential angular velocity 

Angular velocity of the star 

Fluid mass-energy density 

Fluid pressure 

Baryon density 

Fluid 4-velocity in the Einstein/Jordan frame 

Previous quantities evaluated at the star’s center, r = 0 

Stellar radius in the Einstein/Jordan frame 

Asymptotic value of the scalar field 

Gravitational mass of the star 

Scalar charge of the star 

Baryonic mass of the star 

Equation of state parameters, see Eq. (3.171 


Table 1. Variables and conventions used in this paper. Quantities above the horizontal line define 
the theory; quantities below the horizontal line refer to stellar models. 
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is denoted as The energy-momentum tensor of the non-gravitational fields 

is defined by 


2 5Sn^[A^{ip)gp^]'^] 


T = — 

and its trace is given by T = 

The energy conservation equation reads 






T A- T^v) 
A- 


= 0, 


or, more explicitly, 


d log A{<p)^^ 

V -L ^ ^ 


dip^ 


(2.4) 

(2.5) 

( 2 . 6 ) 


Here 


rTi{ip) — 2 5Sip[gpcr\‘p\ 

1 

47rG* 




-V{p)g^ 


fj^y 


(2.7) 


is an effective energy-momentum tensor for the scalar fields, where denotes the 
scalar kinetic and potential contributions to the action (2.11. One may build an energy- 
momentum tensor which is conserved with respect to the Levi-Civita connection of 
the Jordan-frame metric, and whose components are directly related to physically 
observable quantities as 


^ _ 2 SSmlOpa;'^] 


= A-\p)Tp 


( 2 . 8 ) 


2.3. Complexification 

If the target space has a Hermitian structure, then it is useful to write the action in 
terms of holomorphic coordinates and complexified tensors: 


+ {<P, 0)9p.v\ Ik] . 

The complexified field equations are: 


^ - 9^''lab{p, V?) -V{(fi, ip) 


(2.9) 


Rpv =‘i')ab{p,p)^+ ‘^V{p,p)gp^+ -]^Tgp^ , ( 2 . 10 ) 




( 2 . 11 ) 


dip'^ ^ ’ dipf> 

Note that for Kahler manifolds (and in particular one-complex-dimensional manifolds) 


=0, 7V(7>,<^) = , 

so in this particular case the scalar field equations would simplify considerably. 


( 2 . 12 ) 
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2.^. A two-real-scalar model with maximally symmetric target space 

The simplest extension of a ST theory with a single real scalar field is the case of two 
real scalar fields. We will mostly focus on this model to illustrate the basic features 
of the new phenomenology arising in TMS theories relative to the case of a single 
real scalar. If the target space is assumed to be maximally symmetric, then there are 
three possibilities for its geometry: fiat, spherical, or hyperbolic. In the fiat case, the 
target space may be trivially identified with the complex plane C. In the spherical 
case, the target space may be conformally mapped to the one-point-compactification 
C of the complex plane C by means of stereographic projection. In the case of a 
hyperboloid of two sheets, the target space may be conformally mapped to C\S'^, also 
by means of stereographic projection (we shall neglect the case of a hyperboloid of 


one sheet); see Appendix A for details. Using the complex formulation discussed in 
Section |2.3| we work with a single complex scalar rather than two real scalars, for 


which the action (2.91 reduces to 


S = 


dTrG, 


J 


R 

_ - (f) 


+ <S'in[A ((/?, (p)gfj,u', Ik], 


(2.13) 


and the field equations are 




14) 


51og7((/7, (p) 


-47rG*7 ^ ——T. 




dp 


(2.15) 


Hereafter we assume that the potential vanishes, i.e. V{p, p) = 0, and that the target 
space is maximally symmetric. Therefore, upon stereographic projection and field 


redefinition (see Appendix A I, the target-space metric can be written as 

1 / 


7(^,^)=2(l + g) , 


(2.16) 


where r is the radius of curvature of the target-space geometry: for a spherical 
geometry we have > 0, for a hyperbolic geometry < 0, and in the limit r —>■ oo 
the geometry is fiat. 

With the above choices, the field equations become 


= 2 


U\p = 


(l + g) d(^pd,)p + SttG* - -Tg^,J , (2.17) 

g'^'^d^pd^p - 47rG* (^1 -k k { p , p) T, (2.18) 


2p 


pp -\- At^ 


where we introduced 


/ _ o G , W\ d\og A{p,p) 


( 2 , 19 ) 


the so-called scalar-matter coupling function. 
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The function (p), whose derivative enters into the field equations, determines 
the scalar-matter coupling through Eq. (2.19). Without loss of generality we assume 
that far away from the source the field vanishes, i.e. that the asymptotic value of the 
scalar field is ipoo = 0. We then expand the function log A in a series about tp = 0: 


\ogA{ip,ip) = a if + a*if + + • ■ • , (2.20) 

where /3o is real, while a* and /3J are in general complex number^ Although the 
five real parameters Re[Q;*], Im[a*], ,5o, Re[/3J], Im[/3j['] are defined in terms of a specific 
target-space coordinate system, the four real quantities (|a*|, /3o, |/3*|, arga* —^ arg/3J) 
may be expressed solely in terms of target-space scalar quantities, and thus have an 
invariant geometric meaning The remaining real parameter is an unmeasurable 
overall complex phase. 

To make this explicit, redefine = /3ie'®, where 0 is chosen such that jSi is 
real. Then, after defining a* = and a new field ip = the field equations 

become 


= 2 ( 1 + ^ ) + 8ttG. [ 


I 5 


apj = 


24 ! 


ipip -f 4x^ 


tp4’\ 


- 47rG* ( 1 + ^ ) "0) T , 


where the function k is defined in Eq. (2.19) and 
log A{ijj, ip) = aip + ofip + 


( 2 . 21 ) 

( 2 . 22 ) 


(2.23) 


Therefore, any solution of the original theory (formulated with respect to ip and 
complex coupling coefficients a* and PI) can be obtained from a theory where we 
consider the field ijj, a real-valued /3i and a generically complex a. The solution for 


the theory corresponding to the conformal factor (2.20) is then given by a simple 
rotation, p = ip exp (—i0/ 2). 

The model just described represents the simplest, yet quite comprehensive, 
generalization of the model of single ST theory investigated originally in Ref. . 
Note that the quantity \ap = aa = Re[a]^ -I- Im[Q;]^ is strongly constrained by 


observations (cf. [Appendix B I, similarly to the single-scalar case. When a = 0, the 
conformal coupling reduces to 


log A{xp, tp) = ^Pofpip A^Pi'ip'^ ■ 




(2.24) 


where we neglected higher-order terms in the scalar field. However, in TMS theories 
a is a complex quantity and its argument, arga, is completely unconstrained in the 
weak-field regime. In Section |3.2.2| we will show that compact stars in theories with 
a = 0 and a 0 are rather different. 


* At the onset of spontaneous scalarization |(^| <C 1, and we can always expand the conformal factor 
as in Eq. ( |2.20[ |. For scalarized solutions the field amplitude may be large, the higher-order terms 
in the expansion may not be negligible, and the expansion ( |2.20[ | should be considered as an ansatz 
for the conformal factor. For a general functional form of the conformal factor, the series expansion 
used here (and in Ref. [23| ) can only provide a qualitative description of the scalarized solution, 
t The eigenvalues of the quadratic form in (2.20i, given by /Sq di are target-space scalars. The 
phase difference arg a.* — ^ arg 0* arises when this quadratic form is contracted with a* , see Eq. (B.6 i. 
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The field equations can be also written in terms of two real scalars. For this 
purpose, let us split the field ip into real and imaginary parts: ip = Re['!/i] -|- ilm[i/;]. 
Then the conformal factor (2.241, again in the a = d = 0 case, reads: 


\ogA{ip,ip) = ^[{f3o +Pi)Re[ip]'^+ {Po-Pi)^T^bP?] ■ (2.25) 

The structure of this TMS theory is ultimately determined by three real parameters: 
pio+ /3i, Po — Pi and the target-space curvature defined by . When a ^ 0, two further 
parameters (|a| and arga) are necessary to define the theory. 


3. Stellar structure in tensor-multi-scalar theories 


In this section we consider the structure of relativistic stars in the context of the TMS 
theory introduced in Section [2.4[ We first derive the equations of structure for a slowly 
rotating star in the Hartle-Thorne formalism (4^[^ (Section |3.1[ ), then we integrate 
these equations and discuss some properties of scalarized solutions in increasingly 
complex scenarios (Section 3.21. 


3.1. Equations of hydrostatic equilibrium 

We describe a stationary, axisymmetric star, composed by a perfect fluid, slowly 
rotating with angular velocity 12 , using coordinates = (t, r, 0 , (p) and the line element 

guvdx^dx" = ^ -I- r’^{d9^ + sin^ Odp^) 

1 - 2g{r) 

+ 2 [u}{r) — n]r^ sin^ 9dtd(p. (3.1) 

where we neglect terms of order ^ 12 ^ and higher in the metric and in the 
hydrodynamical quantities. The variable /i(r) is related to the more familiar mass 
function m(r) by g = m/r. The energy-momentum tensor of the perfect fluid takes 
the usual form 

= A^{iP, fp) [{p + P)u>^u'' + Pg'^"'] , (3.2) 

where p, P, and = A~^(ip,ip)u^ are the mass-energy density, pressure, and four- 
velocity of the fluid, respectively, and 

0, 0, 12). 


(3.3) 
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With these choices, the field equations (2.2)-(2.3) reduce to a system of coupled 
ordinary differential equations, namely 


1 


= o(l - 2 m)?' 1 + 


'ily'ij; 

4r^ 


'4y)Tp , 

-1 


P' = -{pc^ + p) J ^ + i tpw + K{ip, 


2p 


r(l - 2 m) 






, , 8'!rGi,A'^{'ip,'ip)rP 


V ' V ' + 4 r 2 ?'(1 ~ 2 m ) 


-f 


47rG*y4'‘(-!/','!/') 
1 - 2m 


(m-3P) 


(-i) 


1 - 2m 


V') + rip^p - P) 


OJ = 


, T , yi’ 

’'11 + 55 


'ip'-ij;' - 

r 


, ^ 47rG*54"^(^/’,-i/’)r(pc^-I-P) 
1 - 2m 




(3.4) 

(3.5) 

(3.6) 

(3.7) 


uj' + -uj \ , 


(3.8) 


■where primes denote derivatives with respect to the radial coordinate = r. As 
usual, the system is closed by specifying a barotropic equation of state P = P{p). 

For the purpose of numerical integration, it is useful to work out series expansions 
of the functions p,{r), i'{r), P{r), ip{r), and a;(r) about r = 0: 


M(?’) = ^4Por^ + 0(r"), 

’^{r) = vo + + 3Po)?'^ + C>(?'^), 

P{r) = Po-^Ao(po + Po) [Po + 3Po + (po - 3Po)] r'^ + 0{r ‘^), 

’fpir) = -00 + ^^^A;)ko(Mo - 3Po) (^1 + . 

uj{r) = Wo -k *^^a;oAo(po + To)?"^ + 0{r^). 

5 


(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 


Here the subscript 0 denotes evaluation at the stellar center r = 0. 

When ip and ip are constant and A{ijj, ^/i) = 1 the field equations reduce to the 
standard Tolman-Oppenheimer-Volkoff equations in GR. Indeed, GR solutions are 
part of the solution spectrum of TMS theories with a = 0 in Eq. (2.231, as in the 
usual single-scalar case (T^[^ . On the other hand, when a ^ 0 the scalar field is 
forced to have a non-trivial profile in the presence of matter (T ^0). 

Furthermore, even when a = 0, other solutions characterized by a non-trivial 
profile for the scalar fields can co-exist with the GR solutions; in Appendix G we 


give a simple interpretation of these “spontaneously scalarized” solutions in terms of 
a tachyonic instability of relativistic GR solutions. Besides their gravitational mass 
M and radius R, scalarized neutron stars are characterized by their scalar charge Q, 
which is generally a complex number for the complex scalar field ip discussed her^ 


* Note that Q is defined in terms of a specific target-space coordinate system. Observable quantities 
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The gravitational mass M and the scalar charge Q of stellar models in TMS 
theories can be computed by integration in the vacuum exterior region, where 
P = p = 0. We integrate the structure equations outwards from r = 0 up to a 
point r = R such that P{R) = 0, which determines the stellar surface in the Einstein 
frame. The areal radius R in the Jordan frame can be obtained rescaling R by the 
conformal factor which depends on the value of the scalar field ijj and its 

complex conjugate ip tA r = R. The values of the mass function m, of the scalar 
field Ip and of its derivative ip' at r = R are used as initial conditions to integrate the 
structure equations in the exterior region. In practice, the integration is terminated 
at some finite but large grid point where r = i?co- From the values of m, ip and ip' at 
Roo we can determine M and Q by solving the system of equations 




2 r 


2r2 


Hr) + ^ + 


(3.14) 

(3.15) 

(3.16) 


where ipao is the constant background value of the scalar field. Therefore we have 
three equations to solve for three unknowns: M, Q and ipoo- Th® physical solution 
of interest is the one corresponding to the particular central value of the scalar field 
ipQ (which can be found e.g. by a shooting method) such that the background field 
vanishes, i.e. ipao = 0 . 

As an alternative integration technique, we have also implemented a compactified 
coordinate grid in the vacuum exterior introducing a variable y = 1 /r, which results 
in a regular set of differential equations that is readily integrated to spatial infinity 
at y = 0. The gravitational mass M, scalar charge Q and asymptotic scalar field 
magnitude are then directly obtained from m{y = 0 ), ip{y = 0 ) and dyippy = 0 ), 
and numerical shooting provides a fast-converging algorithm to enforce the boundary 
condition limj,_,,o ip = 0. The two independent integrators yield bulk stellar properties 
that agree within ^ 1 % or better. 


3.2. Numerical integration and results 

In this section we discuss the result of the numerical integration of the hydrostatic 
equilibrium equations in the TMS theory of Section |2.4| Our main interest is to 
understand scalarization in this model, and for simplicity in this paper we will focus 
on static stars. We will use the polytropic equation of state labeled “EOSl” in Ref. j^, 
for which the pressure P and the energy density p are given as functions of the baryonic 
density ns by 

D , ^noiTiB fnsY 

P = Knoms \ — , p = nBmB+K -- — , (4.17) 

\nQ J 7-lV«o/ 

where ng = 0.1 fm“^, itib = 1.66 x 10“^^ kg, K = 0.0195 and 7 = 2.34. Therefore, 
the function p = p{P) can be constructed parametrically by varying tib- 


must be invariant under changes of coordinates, and can thus only depend on target-space scalars such 
as IQp, or \Qa — Qb\^ (for ^ binary system with bodies A and B), or more complicated contractions. 
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3.2.1. The O{2)-symmetric theory In the absence of a scalar potential, the 
gravitational part of the action (2.11 is invariant under the target-space isometry group 


G. For our simple two-real-scalar model with maximally symmetric target space, G is 
the orthogonal group 0(3) in the case of spherical geometry, the indefinite orthogonal 
group 0(2,1) in the case of hyperbolic geometry, and the inhomogeneous orthogonal 
group IO{2) = X 0(2) in the case of flat geometry. 

When scalar-matter couplings are introduced, the action is no longer invariant 
under all of G, but only under some subgroup H < G. As a first example, let us 
consider the particular case in which /3i = a = 0. In this case the conformal factor 


A{'tp,ip) given in Eq. (2.241 reduces to 


A{ip,ip) = exp ( -Mill 


(3.18) 


where again we have neglected higher-order terms in the scalar held. This equation 
is obviously invariant under rotations in the complex plane ("0 —>■ V'e^^) and complex 
conjugation (ip —>■ 'ip). Therefore, H = 0(2). Note that the boundary condition 
ipao = 0 is iJ-invariant. We refer to this special case as the 0(2)-symmetric TMS 
theory. In this theory, a GR stellar configuration with tp = 0 is always a solution that 
is 0(2)-invariant. 

We now construct scalarized solutions, which spontaneously break the 0(2) 
symmetry. They depend on the two real parameters (/3o and r^) of this theory, as 
well as the central baryon density n^. The 0(2)-symmetric character of the scalarized 
solution space is exhibited in Fig. |3.1| where we show that, for given values of r and ub, 
there exists an infinite number of scalarized solutions characterized by a different value 
of the complex held -ipo at the center of the star. The different values of the scalar held 
are related by a phase rotation, and the masses and radii of neutron star models along 
each of the circles shown in Fig. |3.1| are identical. The target-space curvature r has the 
effect of suppressing (r^ < 0) or increasing (r^ > 0) the value of |'0o|j and consequently 
of the scalar charge Q. Therefore a spherical target space (r^ > 0) produces stronger 
scalarization effects in the mass-radius relations with respect to the case of a flat 
target-space metric, as illustrated in Fig. 3.2 On the other hand, a hyperbolic target 
space (r^ < 0) tends to reduce the effects of spontaneous scalarization. This can be 
intuitively, if not rigorously, understood by a glance at Eqs. (2.18) and (2.19): the 


curvature term plays the role of an “effective (field-dependent) gravitational constant” 
which is either larger or smaller than the “bare” gravitational constant depending on 
whether > 0 or < 0. In both cases, as r —>■ oo the solution reduces (modulo a 
trivial phase rotation) to that of a ST theory with a single real scalar held ip and scalar- 
matter coupling A{ip) = exp (i/3ot/’^) ■ We remark that due to the 0(2) symmetry, 
all solutions of this theory are equivalent to solutions with Im[^/>] = 0; as discussed in 
Section |3.2.2| below, these are effectively - modulo a held redefinition - solutions of a 
single-scalar theory. 

Finally, in Fig. |3.3| we illustrate the radial profiles of the mass function m, metric 
potential mass-energy density p and scalar held ip for scalarized stellar models with 
fixed baryonic mass Mb = 1.70 Mq in theories with /3o = —5.0 and = ±1/4. 


3.2.2. The full TMS theory We now turn our attention to the existence of scalarized 
stellar models in the theory defined by Eq. (2.23), which depends on three real 


parameters (/3o, /3i and r^) and the complex constant a. When 0 = 0 and /3i ^ 0, 
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Figure 3.1. Spontaneous scalarization in a TMS theory with 0(2) symmetry. The value i/jq of 
the scalar field at the center of the star for scalarized solutions in the 0(2)-symmetric theory with 
00 = —5.0 and central baryon density ns = 10.4nnuc, where the nuclear density is Tinuc = lO^^m-3. 
Left panel: spherical target space with > 0. Right panel: hyperbolic target space with < 0. In 
both panels the origin corresponds to the neutron star solution in GR. 




Figure 3.2. Stellar properties in the O(2)-symmetric theory. Left panel. The mass-radius relation 
for different values of t and 0o = —5.0. Right panel: Central value of the magnitude of the scalar field 
as a function of the stellar compactness G*M/(.Rc^). Here R is the areal Jordan-frame radius 
of the star. The onset of scalarization does not depend on the value of t. 


this theory is invariant under the symmetry group Z 2 x Z 2 generated by conjugation 
{ijj —>■ '0) and inversion {ip —>■ —ip). Introduction of a S M partially breaks this 
symmetry down to Z 2 , consisting of conjugation only, whereas introduction of a G C\K 
fully breaks this symmetry. 

An interesting question is whether there exists a region of the parameter space 
of this theory in which both fields scalariz^ We first searched for such “biscalarized” 


* This question is not invariant under field redefinitions. More precisely, we ask whether there exists a 
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r [km] r [km] 


Figure 3.3. Radial profiles. Different panels show the profiles of the mass function m, metric 
potential u, total energy density p and complex scalar field -0, in units of c = G = Mq = 1. The 
profiles correspond to a scalarized star in the 0(2)-symmetric theory with = —5.0, /3i = a = 0, 
and fixed baryon mass Mb = 1.70 Mq. The target space curvature is either r = 0.5 (spherical) 
or r = 0.5i (hyperboloidal). In the spherical case, the scalarized solution has a gravitational mass 
M = 1.54 Mq, Jordan-frame areal radius R = 13.0 km, total scalar charge Q = 0.553 Mq and central 
scalar magnitude ji/jol = 0.154. In the hyperbolic case, these quantities are M = 1.54 Mq, R = 13.0 
km, Q = 0.393 Mq and l-^ol = 0.110. For comparison, the GR solution with the same baryonic mass 
has M = 1.54 Mq and R = R = 13.2 km. 


solutions in the Z 2 x Z 2 theory with a = 0, considering a wide range of the (/3o, /3i, r) 
space, but we could not find any. However, the situation is dramatically different 
when a ^ 0. Crucially, |a| has to be small enough to satisfy the observational bounds 
summarized in Appendix B and in particular Eq. (B.3), but arga is completely 


unconstrained by weak-field observations. Our numerical findings are in agreement 
with an approximate analytical model which will be presented elsewhere |51| . In the 
following we discuss the cases a = 0 and a 0 separately. 


Case a = 0.' breaking the 0(2) symmetry down to Z 2 x Z 2 When a = 0 but Pi yf 0, 
we found solutions where only either the real or the imaginary part of the scalar 
field has a non-trivial profile. Therefore, in this case the circle shown in Fig. |3.1 
for the 0(2)-symmetric theory collapses down to four discrete points on the real and 


imaginary axes (cf. Fig. 3.41. 


doubly scalarized solution which can not be described by an effective single-real-scalar theory. 
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Figure 3.4. Symmetry breaking of the space of solutions. When 01 ^ 0, the 0(2)-symmetric 
solution-space analyzed in the previous section (cf. Fig. |3.1| | collapses down to a {Z 2 x Z 2 )-symmetric 
solution-space. This property of the theory is here illustrated for stellar models with the same equation 
of state and central energy density as in Fig. /So = —5.0 and t = 5.0. 


In [Appendix C we perform a linear analysis of the field equations, deriving the 
conditions for scalarization to occur. From Eqs. (C.2|, ( |C.3 1 and (C.5) we expect that 
scalarized models exist if Po + Pi ^ —4.35 when Rep/^] ^ 0, or /3o — /3i < —4.35 when 
instead Im[’!/'] ^ 0. We have checked this expectation by calculating models for the 
parameter sets (i) 1/x — {), Pi —{) and (ii) 1/r = 2, = 0. For each of these cases, we 

have varied the central density from 10“® km“^ to 0.0015 km“^ in steps of 10“^ km“^. 
We applied our shooting algorithm for a scalar field amplitude |^(r = 0)| S [0,1] in 
steps of 0.1, choosing discrete values of the complex phase 0 = 0, 7r/2 , tt, 37 r/2, and 
varying /3o G [—20,3] in steps of 0.01. For all values of the central density and /3o, 
the shooting method identifies one GR solution model with vanishing scalar charge. 
For sufficiently negative /3o, we additionally identify scalarized models. Among these 
models we then identify for a given value of Pd the scalarized model with the lowest 
baryon mass, and thus generate a scalarization plot analogous to Fig. 2 in for ST 
theory with a single scalar field. The result is shown in Fig. |3.5| The small difference 
between the curves for different curvature radius t likely arises from the small but finite 
amplitude of the scalar field appearing in the lowest-mass scalarized binaries, which 
is a byproduct of finite discretization in the mass parameter space. In the continuum 
limit of infinitesimal amplitudes of the scalar field in scalarized models, we expect this 
difference to disappear completely and the dotted and dashed curves to overlap. This 
is indeed supported by an analytic calculation]^ These results confirm the prediction 
of Eq. (C.51 and agree (qualitatively and quantitatively) with the single-scalar case 
shown in Fig. 2 of [46]. 


* This calculation uses Riemann-normal coordinates at (poo in target space, and finds that target- 
space-curvature terms appear in the fiel d equations at third order in the scalar-amplitude expansion. 
Details will be published elsewhere |51|. 
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Figure 3.5. Minimal baryonic mass of scalarized models. The baryon mass of scalarized solutions at 
the onset of scalarization as a function of (i) + 0i for models where Re['0] is non-zero (left panel) 

and (ii) 0o — /3i for models where lm['0] is non-zero (right panel). Each panel contains 6 curves, 
corresponding to the three values of /3i at fixed r = oo (dashed curves) and the same three values of 
I3i at X = 1/2 (dotted curves). The three dashed curves and the three dotted curves, respectively, are 
indistinguishable in the plot and the two families of dashed and dotted curves are only distinguishable 
in the inset, where we zoom into a smaller region. In both panels, the vertical long-dashed curve 
denotes the value = —4.35 above which we no longer identify scalarized models, in agreement 

with Eq. ( |C.5| t. From Eq. ( |2.25[ | it is clear that the natural parameters are Pq -hand po — p\ when 
the theory is written in terms of the real and imaginary part of respectively. 


Indeed, in this case the analogy with the single-scalar case can be made more 
formal. Let us consider without loss of generality the subspace of the solution space 
in which the scalar field is real, i.e. Z = Re[^] ^ 0, W = Im['i/;] = 0. The kinetic term 
can be put in the canonical form by a scalar field redefinition, i.e. 




(3.19) 


where the two fields are related by Z = 2rtan (^)j and —Trr < Z < nt. For \Z\ <C r 
we have 


Z = Z- 

Replacing this Taylor 
the parameters PotPi 


Z3 

12r2 


+ 0(Z5), 


z = z + 


Z3 

12r2 


0(Z5). 


(3.20) 


expansion in 
remain the 


the 


same. 


conformal factor 
In particular, 


(2.24) we see that 
we obtain A(Z) = 


exp 


(/3o + ^i)^V2 


(plus higher-order terms), i.e., the coupling function coincides 
with that of a single-scalar theory with coupling constant /? = /?o + /?i. Thus, as long 
as jZj <C r, the theory with a = 0 is equivalent to a ST theory with one scalar and 
coupling /? = /3 q -I- /3, (or /? = /3o — /3i, in which case only W = Im[i/)] scalarizes). 
Clearly, this proof also includes the limit r —>■ oo, where the solutions reduce exactly 
to those of a single-scalar theory with the identification /3 = /3o + /3i. 

When the condition |Z| ^ r is not fulfilled, the theory is still equivalent to a ST 
theory with one scalar field, but the form of the conformal factor A changes. These 


theories only differ by higher-order terms in the series expansion (2.20), (2.22), which 
are negligible at the onset of the scalarization. 


In Fig. |3.6| we show the mass-radius relation of scalarized neutron star solutions 
in the non-0(2) symmetric theory for different values of r and /3o + /3i. When the 
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Figure 3.6. Mass-radius relations in the full TMS theory. Analogous to the left panel of Fig. |3.2| 
for three values of the curvature radius of the target metric (r = oo, r = 1 and t = 0.5), /3o + /3i = —5 
(left panel) and /3o -h /9i = —10 (right panel). Here we only consider models where Re[i/i] 0. 
The gravitational mass M is shown as a function of the Jordan-frame radius R. For comparison, 
we include in both panels the GR curve. Note the different axis ranges in the two panels. When 
r —>■ oo, the theory reduces to a ST theory with one scalar and effective coupling /3 = /Jq + /3i, and 
the observational constraint /3o + /3l ^ —4.5 is in place [24| . However, such lower bound might be 
less stringent when r is finite. 


coupling is large, we observe that the solutions can differ dramatically from their GR 
counterpart. 


Case a ^ 0; multi-scalarization When ct yf 0, GR configurations are not solutions 
of the field equations. In particular, a constant (or vanishing) scalar field does not 


satisfy Eq. (2.181 when T ^ 0. Therefore it is not surprising that when a 7^ 0 we 
can find solutions with two non-trivial scalar profiles even when /Ig = /3i = 0. A 
more interesting question is whether there are stellar configurations in which both 
scalar fields have a large amplitude. As we have seen, these “biscalarized” solutions 
are absent in the a = 0 case. Here we present preliminary results that demonstrate 
the existence of interesting biscalarized solutions as long as a 7 ^ 0 . 

For concreteness we set |a| = 10“^: such a small value of |a| satisfies the 


experimental bounds discussed in Appendix B (but we have also studied the case 
where |a| = 10 “^, with qualitatively similar results). 

In this preliminary study we vary arga in the range [0, 27r] in steps of tt/6 and 
we set 1/r = 0 (i.e., we consider a fiat target space). A finite target-space curvature r 
does not change the picture qualitatively; a more detailed analysis will be presented 
elsewhere 51 . Our search yields several models with non-zero scalar field, as shown 
in Figs. and |3.8| where dots denote the real and imaginary parts of the central 
value of the scalar field ipQ for which solutions were found. 

For the time being, we wish to remark two very interesting (and perhaps 
unexpected) features of these biscalarized solutions: 


1) Figure 3.7 shows that the solutions are at least approximately 0(2) symmetric 
when Pi ~ |a|, and the 0 ( 2 ) symmetry is broken (the solution “circles” turn 
into “crosses”) when Pi ^ |q;|. The cross-like shape of the scalarized solutions in 
the Re[V’o]j Im[V'o] plane collapses towards a set of solutions on the vertical line 
Re['!/)o] = 0 for the larger values of Pi (bottom panels in Fig. 3.71. This behavior 
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Figure 3.7. Scalar field amplitudes in the full TMS theory - I, Scalar field amplitude at the stellar 
center i/jq for stellar models with fio = —5, |a| = 0.001 and fixed baryon mass Mb = 1.8 Mq. The 
different panels show the solutions found for different values of fii as indicated in each panel. In each 
case, we vary the phase of a from 0 to 27r in steps of tt/O. In contrast to the a = 0 case in Fig. |3.4[ 
the breaking of the 0(2) symmetry occurs gradually as /3i is increased away from 0. 


can be interpreted as an approximation to the spontaneous scalarization for the 
case a = 0 displayed in Fig. |3.5[ and discussed in the text around Eqs. (3.191 


and (3.201. There we observed that spontaneous scalarization of Re[t/>] occurs (in 
analogy with the single-field case) if /3o+ Pi —4.35, and scalarized models with 
a large imaginary part Im['(/;] exist if Po — Pi % —4.35. The biscalarized models 
in Fig. 3.5 have been calculated for fixed Pq = —5. For Pi > 0.65 we therefore 
enter the regime where Po + Pi^ —4.35, and we no longer expect to find models 
with strongly scalarized Re)!/)]. The condition Pq — Pi < —4.35 for scalarization of 
Im[^/>], however, remains satisfied, so that scalarized models should cluster close 
to the Re['0o]^a'Xis. This is indeed observed in the bottom panels of Fig. 3.7 
Note that in this case the condition Pi > 0.65 ^ |a| = 10“^ is satisfied, in close 
correspondence to the case a = 0 of Fig. |3.5[ 


2) Figure 3.8 (which is a “zoom-in” on the top-left panel of Figure 3.7) indicates 
additional fine structure in the space of solutions, with at least three different 
families of scalarized solutions having remarkably different values of the scalar 
field (and therefore of the scalar charge). 
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Figure 3.8. Scalar field amplitudes in the full TMS theory - II. The data of the upper left panel of 
Fig. |3.7| are shown on different scales to resolve the fine structure of the solutions in the i/jq plane. In 
each panel the vertical extent is equal to the horizontal. 


When r —>■ (X), binary pulsar observations in the single-scalar case would impose a 

More work is required to clarify whether a 


constraint equivalent to /3o+/?i ^ —4.5 24 


similar constraint is in place also for multiple scalars. Preliminary calculations indicate 
that the target-space curvature should affect the energy flux from compact binaries at 
high post-Newtonian order. However, it is unclear whether the formalism of for 
describing orbital binary dynamics is applicable to the theory studied in this paper, due 
to the discontinuity at a = 0. Furthermore, for multiple scalars, it is possible that some 
combination of /3o and /3i other than their sum may be constrained by binary pulsar 
observations. A detailed answer to this question requires two theoretical developments 
that are currently missing: (1) the investigation of stellar structure in generic TMS 
theories to understand the effect of the theory parameters on stellar properties, and 
(2) an implementation of these stellar structure calculations in flux formulas similar 
to those derived in 13 (or generalizations thereof). These are important tasks that 


should be addressed in future work. In Fig. |3.6| we adopt an agnostic point of view 
and use large values of /3o + Pi in order to illustrate the effect of scalarization in TMS 
theory in some extreme cases. The phenomenological implications and the stability of 
biscalarized stellar models will be discussed elsewhere (^. 


4. 3+1 formulation of the field equations for numerical relativity 

Studies of the strong-field dynamics of compact stars and black holes, whether isolated 
or in binary systems, require the fully non-linear theory without any symmetry 
assumptions. Such studies can now be carried out using numerical relativity 
techniques, and they have already led to new insights into the behavior of ST theories. 
For example, numerical simulations of neutron star binaries in single-scalar theories 
have identified a new phenomenon (“dynamical scalarization”) occurring in the late 
stages of the inspiral, just before merger 27 52 . Similarly, studies of equilibrium 
sequences of neutron star binaries have shown that dynamical scalarization can lead 
to a reduction of the number of gravitational-wave cycles with respect to GR (^ . 
Scalar radiation has also been identified in the inspiral of black hole binaries when the 
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binary is embedded in a non-trivial scalar background gradient |44[|45| . All numerical 
studies rely on a formulation of the theory suitable for numerical implementation, 
which is most commonly achieved in terms of a space-time (or 3+1) decomposition 
of the field equations. Here we present the 3+1 formulation of the field equations for 
general multi-scalar theories. The work presented in this section is a prerequisite for 
future numerical studies of multi-scalar theories, and we hope that it will motivate 
other researchers to investigate this interesting, unexplored topic. 


4 ..I. Action, stress-energy tensor and field equations 


We consider a multi-scalar theory described by the action (2.11: 
1 


S = 


AttG, 


J d*x^/^ 




(4.1) 

For computational purposes it is useful to consider the scalar fields themselves as 
ordinary matter, described by the stress-energy tensor defined in Eq. (2.7 

= 2jab , (4.2) 

while T^i/, defined in Eq. (2.4|, is associated to the fields dt (for instance, the fluid 
composing a neutron star). The total stress-energy tensor, then, is + Tjff). This 
allows us to use the 3 + 1 decomposition of Einstein’s equations given in |53] , with the 
replacement 


+ Since 

[R - 2-iABg°'^- AV) 


- 87rG*r^i. = 0 , 


(4.3) 


Einstein’s equations have the form — {l/2)g^,jR = 8nGi,(T^i, + The trace 

of Eq. (4.31 yields 


-R + 2'fABg'^d^ip^dpp^ + 8E - SttG^T = 0, 


(4.4) 


where T = and therefore we have 


Rf.. - 2jABd^p'^d,,p^ - 2g^,V - 8^G* ( ) = 0 . (4.5) 


By varying the action (4.11 with respect to one gets the scalar field equation 
5S 


47rG, 


7 




=Up^ + 






(4.6) 


where 'Jcd Christoffel symbols on the target space. 
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J^.2. 3-1-1 decomposition 

As discussed in Ref. (see also [^), we consider a slicing of the spacetime in a set 
of surfaces S. We introduce the normal to those surfaces and the projector 

, ( 4 . 7 ) 


and write the metric in the form (/r, = 0,..., 3; i, j = 1, 2, 3) 

ds^ = g^j^^dx'^dx'' = —a^dt^ + hij{dx'' + P^dt){dx^ + (3^dt) , (4.8) 


where a, /3* and hij are the lapse, the shift, and the metric on E, respectively. We 
remark that in these coordinates = 0, therefore the 3-dimensional metric coincides 
with the projector (4.7) restricted to the spatial indices. It is also worth noting that 
= 0, and = 3. 


We define the covariant derivative on Yi as Di = hi^Va- Since dt = an -|- j3, 
the Lie derivative with respect to n is = [dt — Cp)la. The extrinsic curvature is 
defined as 




(4.9) 


Its contravariant form is purely spatial, i.e., = 0. The extrinsic curvature satisfies 

the relation Kij = Cnhij /2, so the evolution equation for the metric reads 




Other useful relations are 54 


W'^n^, = -K, = D''{\na), 


(4.10) 


(4.11) 


where we defined K = 

In the same way, we can define the curvature of each of the scalar fields as 
= —£„(^"^/2. Consequently, the evolution equation for the scalar fields reads 


Cn^^ = -2K^, (4.12) 

where we note that 

It will also be useful to decompose the quantity g ^^^as follows: 
g°‘^dap^dp(p^ = (/i“^ - n°‘n^)daP^dpp^ = D^ip^Dip^ - . (4.13) 


4-2.1. Einstein’s equations The 3-1-1 decomposition of Einstein’s equations with 
matter is given e.g. in Eqs. (2.4.6), (2.4.9) and (2.5.6) of [^; in those equations the 
matter terms are expressed in terms of the quantities p = = —h’‘^n''Tf 




and Sij = /ij 


’’h-^Tap. We simply replace in those equations —)■ +tI4C\ where 


the explicit expression of is given in Eq. (4.2), i.e., we replace p 
f ^ s,, ^ S,, + where: 


P 


n(V>) 


87rG*p(‘") = 7AB + 4K^K^] + 2V , 

87rG*j*(‘") = -2^abD’p^{-2K’!) = AjabD^P^K, 


= 2jab 




B 


+ 2h,,K^K^ - -/i,,Av? 


(4.14) 

(4.15) 


-2hijV. (4.16) 
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We also have: 


47rG* 


- 2 S\f 


-2^ABDiip^Djip^ - 2hijV. 


Then Eqs. (2.4.6), (2.4.9) and (2.5.6) of give: 

(3)i? + K‘^ - =167rG*p, 

+ 2jab + 4V, 

- R^K) =8TTG.f + A-tabD^V’^K^ 


and 


(4.17) 


(4.18) 

(4.19) 


C K 


D,Dj(\na) + + KK,, - 2K,kK'^^ + 4^G* [%(5 - p) - 25„] 

2jABDi^"^Djp^ - 2hijV , 

(4.20) 


where and ^^^R are the Ricci tensor and the Ricci scalar of the metric hij, 

respectively. 


4-2.2. Scalar field equation To decompose the scalar equation (4.6 1 , i.e. 




..ab 


dV 


ABd^ogA^ 


dp 


B ' 47rG*7 T 0, 


(4.21) 


we start by considering the first term, U\p^ (the single-scalar case was discussed 
in [^[^). We have: 


(4.22) 


Up^ = V^ig'^PVpP^) 

= V, [{-n'^nP + h’^P)VpP^] = Va[2n‘^K^ + D^p^] 
= 2CnK^ - 2KK^ + D,D^p^ + Dp{\na)DPp^ . 

Then, since T = S — p, the scalar field equation takes the form 

C„K} = Kip - \d,D'P - i&PDAna) 

- (D-Pd,P - 4AyO + 

_n.^r ..AB d^ogA ^^ 

27rG*7 [S p). 


(4.23) 


4-3. 3-1-1 equations for 2-sphere and 2-hyperboloid 

Let us now specialize to scalar fields living in a two-dimensional target space with 
maximal symmetry and positive or negative curvature, i.e. a sphere or hyperboloid. 


For simplicity we also assume a vanishing potential {V = 0). As discussed in Appendix 


[A| the sphere and hyperboloid are both described in stereographic coordinates by the 
metric 

' 1 0 


lAB = F 


0 1 


(4.24) 
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with 


F{Z, W) = 


[(^2 -ttT2)/4-hr2] 


2 • 


(4.25) 


Here (p^ = {Z^W), and is positive (negative) for the sphere (hyperboloid). The 
Christoffel symbols (see Appendix AI are 

1 


lAB = 


-2Z -2W 


_ 

I Iab — 


1 f 2W -2Z 
Z2j^W'^+x‘^ I -2Z -2W 


^2 p{/2 ^2 -2W 2Z 

and Einstein’s equations can be written as 

(3)_r + k‘^_ = 2-fAB + 4K^K^] + WirG^p 

= 2F [D^ZD,Z + D^WDiW + 4(A:| -h K^)] + IGirG^p , 

(4.26) 

- F^K) = ijABD^p^K^ + SnG^f 

= 4F{D^ZKz + D^WKw) + SttGJ* , (4.27) 

LnK,, = - D,Dj{\Aa) + + kK,j - 2 K,kK) 

- 2^ABD^p'^Djp^ + 47rG* [%(5' - p) - 2Sij] 

= - AG, (lna) + + kK,j - 2K,kK) 

- 2F{DiZDjZ + D,WDjW) + 47rG* [hij{S - p) - 25'^] . 


(4.28) 


Finally, the scalar field equations are 


CnKz = KKz - \d,D^Z - iG*ZA(lna) - - dAT^A'^) 

= KKz - \d.D^Z - i AZA(lna) + ^ ^ 

X [{ZDiZD^Z - ZDiWD^W + 2WDiZD^W) 

-4 {ZKl - ZK^ + 2WKzKw)] - 2 ttG^F-^^^^{S - p ), (4.29) 


1 


1 


UKw = KKw - - -D^WD,{\na) - - 4K^K^) 


1 


2 

= KKw - \d,D^W - \D^WD,{\na) + + ^2 + ,2 

X [{-WDiZD^Z + WDiWD^W F2ZDiZD^W) 

-4 {-WKl + WKl, + 2ZKzKw)] - 2 t:G.F-^'^^^{S - p) 


(4.30) 

(4.31) 


where we used the expressions of the Christoffel symbols given in [Appendix A[ 


5. Discussion and conclusions 

In this paper we have barely scratched the surface of the potentially rich 
phenomenology of gravitational theories with multiple scalar fields. Several important 
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issues should be addressed in follow-up work. First of all, it is important to compute 
experimental bounds on the parameters /Sq, |/^i|j r, |Qf| and arga— ^ arg/3i that follow 
from binary pulsar timing [T| |24| . The quadrupole-order scalar and tensor radiation in 
TMS theories was computed in |13] , but it is unclear whether the formalism of 13 is 
applicable to the theories that we have studied due to the discontinuity at |a| =0. In 
any case, drawing exclusion diagrams in the multidimensional phase space of the theory 
would require extensive stellar-structure calculations, that will be presented in future 
work. Preliminary results suggest that the target-space curvature radius r enters the 
gravitational-wave flux (at least formally) at high post-Newtonian order, and therefore 
it is quite likely that r will be poorly constrained by binary pulsars. This opens the 
possibility of interesting new phenomenology in the sensitivity window of advanced 
Earth-based gravitational-wave detectors. Furthermore, it is unclear whether binary 
pulsar observations will constrain /3o, |/3i|, or some combination thereof, and the 
parameter arg a — ^ arg /?i (which according to our preliminary results plays a crucial 
role in “biscalarization”) is presently unconstrained. 

In this work we have mainly presented various formal developments, but also 
some new physical results, which in our opinion are representative of the behavior of 
a wide class of TMS theories: 


(i) In theories with a = 0, GR solutions co-exist with scalarized solutions but 
(besides the case of 0(2)-symmetric theories with Pi = 0) we could not find 
any “biscalarized” solution for any value of Pq. In other words, in this case either 
the real or the imaginary part of the complex scalar field scalarizes but not both, 
and the 0(2) symmetry of the Pi = 0 case is broken. 

(ii) The a yf 0 case is dramatically different. In this case - even when |q;| is small 
enough to be compatible with Solar System constraints - biscalarized solutions 
exist, and their existence depends quite critically on the value of arg a — ^ arg Pi. 
These solutions seem to exist quite generically in the complex-a plane, but their 
properties strongly depend on the values of Po and \Pi\. 


These results were obtained through extensive numerical searches. However, given the 
large dimensionality of the parameter space, we cannot exclude the existence of other 
solutions which were not found in our initial searches. An approximate analytical 
model which supports our results and a more detailed analysis of biscalarization will 
be presented elsewhere [sT] . 

Some obvious extensions of the present results concern the study of rotating 
neutron stars in TMS theory (generalizing (4^1 |56]-|^ ) and of the universal relations 
valid for neutron stars in general relativity which may or may not hold in 

this theory. Another possible extension is to relax the assumption of a vanishing 
potential in the action, i.e. to consider the multi-scalar generalization of massive 
Brans-Dicke theory 61,62, and investigate the effect of the scalar field masses on 
the structure of scalarized neutron stars. The use of more realistic equations of 
state is pivotal in confronting TMS theory predictions on the evolution of binary 
pulsars with observations. Furthermore, we hope that the 3-1-1 split worked out 
in Section ^ will encourage other research groups to perform numerical simulations 
in TMS theories. This may lead to studies of phenomenological interest, such as 
investigations of dynamical multi-scalarization in neutron star binaries and evolutions 
of binary black-hole systems in the presence of non-trivial scalar field backgrounds. 
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Appendix A. Spherical and hyperboloidal target spaces 


In TMS theory [^, the scalar field is an application from the space- 

time manifold At to the target-space manifold T. This target-space manifold is 
Riemannian, and its metric is denoted by The dimensionality of T (i.e. 

the number of scalar fields) is N. Since one-dimensional manifolds are necessarily flat, 
the simplest non-trivial case is N = 2. Furthermore, the simplest two-dimensional 
manifolds are the maximally symmetric ones, i.e. spherical, hyperbolic and flat spaces. 
In these spaces, the curvature radius f > 0 is constant; the Ricci scalar is 77 = 2/r^ 
for spherical space, 77 = —2/f^ for hyperbolic space, and 77 = 0 for flat space. For 
convenience we define r = f, ir for spherical and hyperbolic spaces, respectively, so 
that the Ricci scalar has the form 77 = 2/r^ in both cases. 

Here we derive the expression for the target-space line element = 

2^d(pd(p in terms of the complexified scalar field tp = Z + iW for the spherical and 


hyperbolic cases; the result is Eq. (2.161 in the main text. 


Spherical target space 

The 2-sphere can be defined from its embedding in a three-dimensional Euclidean 
space of coordinates {x, y, z) through the equation 

x"^+ 1/^ + . (A.l) 

It can be parametrized in polar coordinates, defining (p^ = (0, $), where: 

X = r sin 0 cos 4>, y = r sin 0 sin $, z = rcos0. (A.2) 
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The target-space metric in these coordinates is 


lA'B' 



(A.3) 


This frame has undesirable features: in the flat-space limit r —)■ oo the metric diverges, 


and the kinetic term in the action (4.1 \ [where for simplicity we set 1^(0, 4*) = 0] 


S = 




f - {d^ed^e + sin^ 


(A.4) 


diverges as well. Moreover, the polar frame is not the most suitable for numerical 
implementation, because it has coordinate singularities at the boundary of the open 
intervals O< 0 < 7 r, 0<$<27r where the coordinate system is defined. To 
fix this problem we perform a stereographic projection from the north pole of the 
sphere (which is the only point of the manifold not covered by this chart) to the plane 
= {Z, W) tangent to the south pole: 

2 r sin 0 , . , 

Z=- - x = 2r- -—cos 4), (A.5) 

r — z 1 — cos 0 

Or 9111 

W=- - y = 2f ^ ^ sin$. (A. 6 ) 

f-z^ 1-COS0 ^ ^ 


With this projection the equator is mapped to the circle Z‘^ + W‘^ = 4r^; the upper and 
lower hemispheres are mapped to the exterior and interior of this circle, respectively, 
and the north pole is mapped to infinity. Using sin 0/(1 — cos0) = cot(0/2), the 
complex held Lp = Z + iW is written in a more compact form as 


<p = 2 r cot( 0 / 2 )e'^ . 


(A.7) 


In the coordinate frame = (Z, W) the target-space metric is 

t-i 


(1 — cos 0 )^ 

lAB = " -- 'OAB = 


[(Z2-tlU2)/4-fr2]- 


t 5 ab ] 


(A. 8 ) 


note that {Z'^ + f"^ = 2r^/(l — cos0). In terms of the complex held p, 

SABdp^dtp^ = d(pd(p, therefore "fABdip^dp^ = 2^dpd(p (note that we denote 7 = 7 ^ 5 , 
because a,b can only take the value 1 ), and 



(A.9) 


Hyperbolic target space 

The two-dimensional hyperbolic space of two sheets can also be defined from its 
embedding into with coordinates (x,y,z) through the equation 

-x^ - y"^ + z"^ = r'^ . (A.10) 


It can be parametrized in terms of = (0, 4>) as 
X = f sinh 0 cos $ , y = f sinh 0 sin $ , 


z = ±f cosh 0 . 


(A.Il) 
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The target-space metric in these coordinates is 


lA'B' 


r2 0 \ 

0 sinh^ 0 j 


(A. 12) 


As in the case of spherical space (see above) the metric diverges when f —>■ oo, and 
the kinetic term in the action diverges as well. Therefore we perform a stereographic 
projection from the point at the top of the lower branch to the plane = {Z, W) 
tangent to the bottom of the upper branch. With this projection, the upper branch is 
mapped to the interior of the circle = 4f^, and the lower branch is mapped 

to the exterior of this circle. The stereographic mapping reads 

(A.13) 
(A. 14) 

for the upper branch, and 


^ 2f „sinh 0 
Z = - -X = 2r ---- cos $ . 


z + r 


cosh 0-1-1 


X,. 2f sinh0 . ^ 

W = - -y = 2r --r sm $ . 


z + r 


cosh 0-1-1 


Z = 
W = 


2r sinh 0 

-z-f^ ^ ^'’cosh0-l 

2f sinh 0 

- y = 2r - 

—z — f cosh 0 — 1 


cos$, 
sin$ 


(A.15) 

(A.16) 


for the lower branch. The complex field p = Z + iW is then p = 2f tanh(0/2)e*^ for 
the upper branch, and p = 2f coth(0/2)e'* for the lower branch. 

In the coordinate frame p^ = {Z, W) the target-space metric is 


lAB 


(1 ± cosh0)^ 
- ~A - °ab 


[-(^2-^ W2)/4-kf2]' 


:5aB 


(A. 17) 


where the upper (lower) sign refers to the upper (lower) branch; note that —{Z^ -\- 
lT^)/4 -I- = 2r^/(l ± cos 0). In terms of the complex field p, then, the target-space 

metric is 2jdpdp with 



(A. 18) 


Field equations for two-dimensional spherical and hyperbolic spaces 


In summary, the expressions (A.91, (A.181 for the target-space metric in the (two- 
dimensional) spherical and hyperbolic cases can be written in the form of Eq. (2.16), 

i.e. 


7=l(l + 


4r2 / 


(A.19) 


where r = r for a spherical space, and r = if for a hyperbolic space. In the coordinate 
frame p^ = {Z, W) the target-space metric is 


lAB = 


[{Z'^ + W^)/4- 


t^ab 


(A.20) 


for both the spherical and hyperbolic space. Therefore, Eqs. (A.191, (A.201 describe 
a spherical space if > 0, an hyperbolic space if < 0. The limit r —)■ oo yields 
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flat space. If r —oo and the scalar held is restricted to real values, one recovers the 
single-scalar case. 

The Christoffel symbols are: 


l^zz = - 

2Z 

7 zz — 

2W 

X'^+Z'^+W^ ’ 

x^+z^+w^ 


2W 

^Z _ 

7 ww — 

2Z 

X'^ + Z'^+W'^ ’ 

x^+z^+w^ 


In terms of the complex held Lp, writing explicitly 
only take the value 1) we get 


ryZ __ 21 V_ 

/ zw “ r2+z2_|_^^2 , 

(A.21) 

__ 2 Z_ 

f ZW ~ x'^ + z'^+w'^ • 

the indices a, b in (which can 


'yVb = x5^iog7 = - 




4r^ -|- (fiip ’ 


'rVb = 7;9v'^ogj = - 




-|- (ftp 


(A.22) 


The Ricci tensor and Ricci scalar of the target space are TZab = ^~^Sab and TZ = 2r“^, 
respectively. 

Replacing the expression of the metric (A. 181 and of the Christoffel 
symbols (A.22I in Eqs. (2.141 and (2.151 with V{p) = 0 we And the held equations 
for a maximally symmetric two-dimensional target space, i.e. Eqs. (2.171 and (2.18). 


Appendix B. Solar System constraints 


The weak-fleld limit of TMS theories has been worked out in . Specializing these 
results to the theory constructed in the body of the text, and rewriting them in 
complex notation, one finds that the gravitational constant measured in a Cavendish 
experiment is given by 


G — G*A^(1 -I- Koodoo) , 


(B.l) 


where the subscript oo denotes evaluation at pco = 0 and we defined the complex 
function K{ip,ip) as in Eq. (2.19). Using Eq. (2.20), one finds that Kqo = 2a*. 

It is straightforward to show that the post-Newtonian parameter yppN reads |13| 

7PPN -1 = - r.T T = 2 ^ (b-2) 

implies the constraint 


>o _ 8|a*p 

l-kitooKoo l-l-4|a*|2 

and therefore the Cassini bound |7 ppn ~ 1| < 2.3 • 10“® 

< 3 - 10 


22 


(B.3) 


on the coupling constants a* and a* appearing in Eq. (2.201. Crucially, the one above 
is a bound on |q;*|, whereas arga* is completely unconstrained in the weak-fleld limit. 


On the other hand, the post-Newtonian parameter /3ppN reads 13 

^ 00 ^ 00 /^oo 


/3ppN — 1 — 


2(1 -)- Avqq/Coo)^ 


f^OO ^OO 

where the real-valued function (p) is defined by 

1 




d 

d(p 


(B.4) 


(B.5) 

























Tensor-multi-scalar theories 


30 


Using the definitions above and Eq. (2.201 we obtain 


/3ppN — 1 — 


a*a*PI + 2a*a*Po -t a*a*P^ 


2|a*p 


(1 -I- 4|a*p)2 


( 1 - 1 - Aa*a*Y 

Po + |/3i|cos(2arga* -arg/3*) 


(B.6) 


Finally, the bound |/3ppN —1| < 1-1 • 10“^ coming from the combination of Cassini 
and Lunar Laser Ranging measurements 63 implies a constraint on some combination 
of the parameters /3o, \Pi\, and arga* — ^arg/3*. However, note that if |a*| —>■ 0 the 
observational constraint I^SppN — 1| < 1.1 • 10“'* is satisfied for any value of /3o, \Pl\ 


and arga* — ^ arg/3*, and therefore these parameters are unconstrained by weak-field 
observations in this limit. 


Appendix C. Linearized field equations and scalarization 


Here we consider the ST theory defined by Eqs. (2.211 and (2.23) with a = 0, which 


admits GR solutions with tp = 0. We will perturb these GR solutions, and linearize 
the field equations in the perturbations. This is valid when the amplitudes of the 
scalar fields are small and consequently the metric back-reaction on the scalar field 
can be neglected. This approximation is well motivated at the onset of scalarization. 

The field equations acquire a particularly simple form when linearized to first 


order in Z = Re['(/'] and W = Im[■)/>]. In this case, the tensor field equations (2.21) 
reduce to 

~ 5 (C.l) 




and therefore the background geometry to 0{Z,W) is described by a GR solution. 


The scalar-field equation (2.22) becomes 


\3Z — —47rG*(/3o -I- Pi)TZ , 
nw =- Pi)TW, 


(C.2) 

(C.3) 


where, in this perturbative expansion, the box operator is evaluated on the GR 
background solution and the trace of the perfect fluid energy-momentum tensor T 
attains its GR value, i.e. T = —(87rG*)“*i? = 3P — p. Note that the equations for 
Z and W decouple in this limit, reducing to the same equation as in the single-scalar 
case, Gup = —47rG*/3T(p, but with effective coupling parameters /3 = /?o + /?! and 
/3 = /?o — /3i, respectively. 

In the case of a single scalar, the term on the right-hand side of the scalar equation 
can be interpreted as an effective mass term (cf. e.g. [^) 

mis =-AttG^PT . (C.4) 


Because in typical configurations T ~ —p < 0, the effective mass squared is negative 
when P < 0. This signals a possible tachyonic instability which is associated with an 
exponentially growing mode and causes the growth of scalar hair in a process known 
as spontaneous scalarization |23| , as discussed in the main text. In the case of static 
compact stars, it turns out that this instability occurs for P < —4.35, the threshold 


value depending only mildly on the equation of state 64 ■ 66 
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The same reasoning can be applied to Eqs. (C.2) and (C.3l. Because the latter 
are completely equivalent to two copies of a single-scalar equation, scalarization is 
expected whenever 


/3o + ^1 < -4.35 or (3o - Pi < -4.35. 


(C.5) 


Note that these conditions were derived assuming that each scalar field acquires a 
non-vanishing expectation value independently and by perturbing a static GR solution. 
In particular, they do not imply that both fields scalarize when both conditions (C.51 
are satisfied. In fact, biscalarization can be investigated in this perturbative framework 
by studying the linear perturbations of (say) the scalar field W on the background of 
a previously scalarized solution where Z has a non-trivial profile. 
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